Saturation field of frustrated chain cuprates: broad regions of predominant interchain 
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An efficient and precise thermodynamic method to extract the interchain coupling (IC) of spatially 
anisotropic 2D or 3D spin- 1/2 systems from their empirical saturation field H s (T — 0) is proposed. 
Using density-matrix renormalization group, hard-core boson, and spin-wave theory we study how 
H s is affected by an antiferromagnetic (AFM) IC between frustrated chains described in the J1-J2- 
spin model with ferromagnetic 1st and AFM 2nd neighbor inchain exchange. A complex 3D-phase 
diagram has been found. For Li2Cu02 and Ca2Y2CusOio,we show that H s is solely determined 
by the IC and predict H s « 61 T for the latter. Using H s ~ 55 T from our high-field pulsed 
measurements one reads out a weak IC for Li2Cu02 close to that from neutron scattering. 
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Since real spin chain systems exhibit besides a signif- 
icant in-chain coupling also an interchain coupling (IC), 
one may ask: in which cases is this relatively weak IC still 
important or even crucial? From the Mermin- Wagner 
theorem its decisive role for the suppression of fluctua- 
tions is well-known. The IC leads to magnetic long-range 
order (LRO) at T = in 2D and at T < T N in 3D 
(see e.g. Ref. 0). Often one is faced with a situation that 
the (large) in-chain couplings are known with reasonable 
precision, e.g. from band-structure calculations, inelastic 
neutron scattering (INS) or susceptibility data 0, but 
precise values for the tiny (nevertheless important) IC 
are lacking. If frustration is absent, the IC can be deter- 
mined quite accurately, e.g. from Tn analyzed by Quan- 
tum Monte Carlo based studies [U and more approxi- 
mately using mean- field theory for a 3D IC @. But how 
to extract from experimental data a small IC for frus- 
trated systems with weakly coupled chains where these 
methods do not work? Here we address such a 2D/3D 
problem for the case of frustrated spin- 1/2 chains with 
ferromagnetic (FM) nearest-neighbor and antiferromag- 
netic (AFM) next-nearest-neighbor exchange described 
by the spin isotropic Ji-J2-model (IM). Nowadays it is 
the standard model for edge-shared chain cuprates (see 
e.g. 0). This 1D-IM attracted special interest |§-|lj| due 
to a rich phase diagram with multipolar (MP) phases de- 
rived from multi-magnon bound states (MBS) in high 
magnetic fields jl4l4l6l | . Additional AFM degrees of free- 
dom enhance the kinetic energy of magnons. Thus, AFM 
IC might disfavor multi-MBS. A precise knowledge of the 
magnitude of the IC is therefore a necessary prerequisite 
to attack the multi-MBS problem, including a possible 
MBS Bose-Einstcin condensation. Hence, knowing it is 



of general interest [17H20| . 

In this context, Li2Cu02 (see Fig. 1) is one of the 
best studied frustrated cuprates and therefore well 
suited to compare theory and experiment. In particular, 



using spin-wave theory (SWT) the main in-chain and IC 
J's were extracted from INS data and a specific AFM 
IC was found crucial for preventing spiral order in the 
3D ground state (GS) 0. If the saturation field H s can 
be determined experimentally, an excellent opportunity 
occurs to check the INS based IC. But for Li 2 Cu02 
H s has not been measured so far. Here we report 
high- field magnetization data to fill this gap. Our paper 
is organized as follows. First we revisit the ID case and 
provide details of the density-matrix renormalization 
group (DMRG) technique involved. Then we report 
results for coupled chains, among them a rather complex 
phase diagram, and compare our findings with our 
experimental _ff s -data for Li2Cu02 and arrive at an 
almost perfect explanation in terms of a predominant IC. 

We apply the DMRG method [5l| with periodic 
boundary conditions (PBC) in all directions. Seemingly, 
this method is much less favorable for D > 1; however, on 
current workstations using highly efficient DMRG codes, 
spin systems with up to about VlO x \/l0 x 50 sites, 
i.e. 10 coupled chains of length L ~ 50, can be studied. 
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FIG. 1: (Color) Left: Crystal structure of Li2Cu02 with two 
CUO2 chains per unit cell along the 6-axis. Right: View along 
the c-axis on the a&-plane. The main in- and interchain cou- 
plings Ji t 2 and Jic, J/ c : arcs and dashed lines, respectively. 
The normalized ICs read /3i = J[ c /\Ji\ and /?2 = Jic/|Ji|. 
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FIG. 2: (Color) Magnetization vs. external field (DMRG 
data) for different n as shown in the inset, for a — J2/IJ1I = 
1/3, P2 = 0.05, Pi = 0, and L = 24 sites in each chain. Lower 
inset: The 3D arrangement of chains used for the DMRG 
study reported here. Upper inset: Magnetization curve inside 
the 3-MBS region for P2 = 0.005 (blue curve, x) compared 
with that for the 1-magnon C-phase for P2 = 0.05 (red curve, 
o). Note the three times larger step for P2 = 0.005. 

Thus, by taking a proper arrangement of the chains 3D 
lattices can be simulated, cf. the inset in Fig. [2j Let us 
now describe how the block states are constructed: in the 
n x L cluster, where n denotes the number of chains and 
L is the chain length. If we regard n sites in the ac-planc 
as a "unit cell", the system can be treated as an effec- 
tive ID chain with L sites (step 1). This enables us to 
use an appropriate ID array for the construction of the 
PBC (see Fig. 1 of Ref. [22h . In the second step the sites 
within each "unit cell" are arranged into numeric order 
as shown in the inset of Fig. [2J This way, the distance 
between most separated interacting sites can be held at 
11 and 23 in the 4- and 8-chain systems, respectively. 
Since the exchange interactions run spatially throughout 
the system, the wave function converges very slowly with 
DMRG sweep but without getting trapped in a 'false' 
GS. We typically kept m ps 1600 — 4000 density-matrix 
eigenstates in the DMRG procedure. About 20 — 30 
sweeps are necessary to obtain the GS energy within 
a convergence of 10 _7 |Ji| for each m value. All cal- 
culated quantities were extrapolated to m — > 00 and 
the maximum error in the GS energy is estimated as 
AE/\J\\ ~ 10 -4 , while the discarded weight in the renor- 
malization is less than 1 x 10~ 6 . For high-spin states 
[Stot ~ (nL — 10)/2)] the GS energy can be obtained with 
an accuracy of AE/\ Ji\ < 10~ 12 by carrying out several 
thousands sweeps even with m w 100 — 800. Then, we 
obtain the reduced saturation field h s = g[i B H s /\Ji\ with 
high accuracy (e.g. 12 digits as compared to exact solu- 
tions available in some cases). The phase assignment in 
the ID-diagrams shown in Refs. 0, 0-O and in 3D (see 
Fig. 3) stems from an analysis of the magnetization curve 
just slightly below h s as shown in Figs. 2 and 6. The sig- 
nature of each region is the height of the magnetization 



TABLE I: Saturation field h B at a = 0.332, p 2 = 3/76, and 
Pi = 0. Jic has been multiplied by 4 for 2-chain systems. 



L single chain 2-chain system 8-chain system 

16 0.0610480058942 0.315789473684 0.315789473684 

48 0.0616910423378 0.315789473684 0.315789473684 

96 0.0616910487270 0.315789473684 

144 0.0616910487247 0.315789473684 



steps AS' Z =1,2,3, ... for dipolar (1-magnon), quadrupo- 
lar (2-MBS), octupolar (3-MBS), and hexadecupolar (4- 
MBS), etc. phases, respectively. For the quadrupolar case 
the DMRG results are confirmed by the exact hard-core 
boson approach (HCBA) [il| . 

We start with the ID problem. Using its frustra- 
tion parameter a=J2/\J\\, the critical point is given by 
the level crossing of a singlet and the highest multiplet 
state at a c = 1/4. Approaching a c from the spiral side 
a > a c , it is clear that if s 1D (a) decreases and vanishes 
at a c > a > 0. The curve -ff 1D (a) is not smooth: it con- 
sists of quasi-linear parts with an infinite number of slope 
jumps at the endpoints of each quasi-linear part. These 
intervals become shorter and shorter when a — > a c (see 
Fig. S3 of Ref. 0). These non-analytic endpoints reflect 
the changes of multi-MBS related low-energy excitations. 
This specific behavior persists also for D — 2, 3 (see Figs. 
4, 5) and is a signature of quantum phase transitions. 
Concerning Li2Cu02 we stress that a ID-approach yields 
for a = 0.332 @|, h s = 0.0616916 (see Table I) where 
g = 2 and J] = 228 K have been used. Thus, the empiri- 
cal value H s = 55.4 T reported below is strongly underes- 
timated by a ID-approach. Hence, let us consider, what 
happens, if an IC of the type shown in Fig. 1 is switched 
on. In general, a complex phase diagram in terms of 
the IC and a has been obtained (see Fig. 3 and Ref. 0) . 
First, with increasing IC above a critical J^'~(a) the MP 
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FIG. 3: (Color) Part of the 3D phase diagram around 
Li2CuQ2 in terms of P2 and a. h s is given by contour lines. 



2 



phase is removed in favor of one of the two incommensu- 
rate (INC) dipolar phases where H s (given exactly by the 
SWT) is still affected by both in-chain and IC although 
the contribution of the former for a < 0.54 is signifi- 
cantly reduced as compared to the ID case (see Figs. 4, 
5). But at stronger IC for a < 1 one reaches a second 
critical point ^ r ' 2 + where the INC phases are suppressed 
in favor of a commensurate (C) phase with FM in-chain 
correlations (see Figs. 3-5). Then H s depends solely on 
the IC: 



0.4 



gH B H B = N ic ( J ic + J/J for J ic > J- 



cr,+ 



(1) 



where iV; c is the number of nearest interchain neigh- 
bors (8 in case of Li2Cu02). For a > 1 there is no 
C-phasc. Now we take J/ c = for the sake of sim- 
plicity. Taking into account it is possible but electronic 
structure calculations suggest J[ c / J- lc < 0.1, only. Then, 
J"' + (a) = (4a - 1)| Ji |/9, if a < 0.57 which is relevant 
for Li2Cu02- For the remaining interval 0.57 < a < 1, 
the INC2-C transition is of 1 st order, see Ref. M. Turning 
to Li2Cu02 and ignoring a small anisotropy [3j, we have 
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FIG. 4: (Color) Upper: Approximate relative in-chain contri- 
bution to the saturation field measured by its deviation from 
Eq. (1) for different a values vs. 02 in units of its critical value 
/?2 r,+ (see Figs. 3, SI, and text). The yellow stripe highlights 
the region of predominant IC. The full line parts of a curve 
(below the symbol o) can be determined within the SWT, 
whereas the branches given by symbols are obtained by the 
DMRG and the HCBA methods. Lower: in-chain and IC 
spin-spin correlation functions vs. applied field, given by the 
filled symbols ■ and •, respectively, as obtained by DMRG. 
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FIG. 5: (Color) The saturation field h s at a = 1/3 vs. 02 from 
DMRG method. Notice the two critical IC values = 
J£'~I\M = 0-0109 and#; r ' + = J"' + /|Ji| = 0.03704 denoted 
by dashed vertical lines. For INC1- and C-phases, see text. 
Inset: the weak non-linearity of h s vs. 02 ■ 



J ic w 9.04K > J"' + (0.332) = 0.0364|Ji| w 8.2 K [23j. 
Since in the C-phase H s depends solely on J lc , the IC can 
be directly read off from its measured value H S =55A T 
yielding Jw. = 9.25 K very close to 9.04 K from zero-field 
INS data [3j. In the INC1 phase J lc dominates H s . There 
above J"~ (0.332) 0.0109|Ji| = 2.49 K only INC 1- 
magnon low-energy excitations exist. Below J^'~ 3- MBS 
are recovered as low-energy excitations. The transition 
from the 3-MBS- to the INCl-phase is 1 st order. The gen- 
eral situation is shown in Fig. 4. Well below Jf c r ' + the 
saturation field significantly depends on a. The yellow 
stripe highlights the region where the IC is predominant 
as addressed in our title. Here this dependence is weak 
and smooth. Subtracting their classical value, the spin- 
spin correlation functions show that the in-chain fluctu- 
ations vanish much faster than the IC ones for H — > H s 
(see Fig. 4). This explains the surprising result of Eq. (jTJ) , 
too. For systems in the FM subcritical in-chain frustra- 
tion regime a < 1/4 the in-chain contribution vanishes 
by definition. Hence, the external field has to overcome 
the AFM IC, only, and can be used to extract the IC ex- 
actly, if the corresponding g- factor is known (see Eq. (1)). 
This case is realized e.g. for Ca2Y 2 Cu50io [25j. It has a 
similar IC geometry as Li2Cu02, but a 2D arrangement 
of chains, i.e. N- lc = 4. This reduced Nic is overcompen- 
satcd by a larger IC J ic + J/ c = 26 K [2^|. As a result 
we predict H s pa 61 T and refine a recent overestimated 
value of 70 T from extrapolated low- field data [25| . 

Pulsed-field magnetization studies have been per- 
formed at the Dresden High Magnetic Field Laboratory 
in fields up to 60 T. The results taken at T = 1.45 K 
for H || &-axis on a Li2Cu02 single crystal from the 
same batch as in the INS-study Q is shown in Fig. [S] 
The data imply a quasi-linear increase of the magnctiza- 
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FIG. 6: (Color) Left: Magnetization M (left axis) and the 
derivative dM/dH (o) from pulsed fields. Red (green dashed) 
line: auxiliary line to make the upturn near H s better visible 
(classical curve (CC): M/M B = H/H s ). Right: Deviations 
from CC: DMRG vs. experiment for different IC /?2 and a H a 
value within the experimental error bars. Inset: Behavior just 
below H s for the two /^-values shown in the main part. 



tion M(H) between 10 and 30 T, i.e. dM/dH = M' ps 
const. Above about 50 T M' increases notably and pro- 
nounced peaks develop at 55. 4± 0.25 T and 55.1 ±0.25 T 
for two pieces of our single crystal. The sharp drop 
of M' towards at higher fields justifies to attribute 
the peaks with H s . The saturation moment amounts to 
M s ps 0.99±0.06Ai B /f-u. using g b = 1.98 ±0.12 in reason- 
able agreement with g b = 2.047 from low-field ESR-data 
at 300 K 0]. In Fig. [6] we compare M/M s from the 
DMRG with the experiment for Li 2 Cu0 2 . The DMRG 
description in this standard plot is rather good, but it is 
inconvenient to extract parameter values. This is due to 
the weak quantum fluctuations in Li 2 Cu02- The function 
f(M) = H/H s — M/M s reflects the remaining quantum 
fluctuations much better (see Fig. 6). Notice the en- 
hanced fluctuations (i.e. smaller M(H)) for the weaker 
IC (see inset). Noteworthy, the 'width' of the DMRG 
curves stems from the finite steps of M{H) due to the 
finite L. The almost straight shape of the M(H) curves 
shown in Figs. 2 and 6 evidences the 3D character of 
Li2Cu02 in accordance with the large local magnetic mo- 
ment ~ 0.93fi B observed in the ordered phase at low T 
at ft. = 0. The measured H s of our single crystal results 
in J ic = 9.25 ± 0.04 K, close to the INS-data. Finally, we 
note that linear relations of two experimentally accessible 
thermodynamic quantities, the Curie- Weiss temperature 
and H s , yield useful constraints for the in-chain J's Q. 

To summarize, the crucial role of AFM IC in frus- 
trated quasi- ID systems, such as Li2Cu02, for their be- 
havior in external fields and, particularly, the strength 
of H s has been demonstrated. Extracting J lc for that 
system from pulscd-ficld data, we arrived at 9.25 K, very 
close to a previous INS study which confirms the validity 



of the adopted spin model. To extract J lc a _ff s -study 
is preferable over an INS, due its smaller error bars, the 
possibility to work with small single crystals and its much 
higher efficiency (concerning both time and costs). The 
large H s from completely different studies discards any 
ID scenario for Li2Cu02, even more, H s itself is (within 
the isotropic model) independent of the in-chain cou- 
plings J\ and Ji. Thus, our results show exactly that 
for a relatively wide interval < a < 1 and collinear 
magnetic order at H =0, H s depends only on the IC, 
irrespectively of its strength. A complete study of the 
entire phase diagram including the INC2-phase which at 
present is rather of academic interest, only, will be pre- 
sented elsewhere. The MP-phascs from ID studies are 
very sensitive to the presence of IC in the 2D/3D counter 
parts. In particular, they can be readily eliminated by a 
weak AFM IC. Instead new incommensurate phases may 
occur. A study of other systems and their H s within the 
approach proposed here is promising and in progress. 
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RI615/16-1), PICS program (contr. CNRS 4767, NASU 
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In the present Supplementary part we provide the reader with the general phase diagram of the 3D model with 
a diagonal antifcrromagnctic (AFM) interchain coupling (IC) J lc that couples each site in a frustrated J1-J2 single 
chain to its next-nearest-neighbor (NNN) sites in the adjacent chains. This type of IC is realized in Li2Cu02, see 
Fig. 1 of the main text. Its rigorous treatment leads for the incommensurate parts INC1 and INC2 of the general 
phase diagram (see Fig. SI) to somewhat tedious analytical expressions for the saturation field H s provided here for 
an interested reader explicitly. Then, in the spirit of our proposed efficient thermodynamic approach, we present 
simple but exact equations where H s together with the Curie- Weiss temperature can be used to get insight into the 
microscopic exchange couplings. Finally, we present DMRG data mentioned in main text for the most realistic case 
of a very weak nearest-neighbor (NN) IC /?i = J[ c /\ J\\ and a predominant NNN IC /?2 = Jic/\Ji\ with parameters 
relevant for Li2Cu02 suggested by band structure calculations [lj. 

PHASE DIAGRAMS AND THE DETERMINATION OF THE SATURATION FIELD 

We start with the presentation of the whole phase diagram see Fig. SI. 

Let us first remind the behavior of the saturation field at the basis line of the general 3D phase diagram, i.e. the 
ID limit of the considered problem, see Fig. S2 and also Refs. 15 and 16 in the main text. At present the higher 
multipolar phases beyond the four-magnon sector near the critical point a c = 1/4 cannot been assigned numerically 
with sufficient precision as a sequence of 5-, 6-, etc. phases as suggested by Krivnov and Dmitriev on general grounds 
The question mark in Fig. S2 reminds the reader of this numerically still somewhat unclear situation. 

The hard-core boson approach provides an exact solution (see e.g. Rcf. 11 in the main text) for the 2-magnon 
(nematic, quadrupolar (QP)) phase. In the ID case we obtain in the QP phase 

/i B = 2a-l + 0.5/(a+l), (SI) 

valid above a critical point a 3 w 0.3676776, i.e. for a > CK3. Below a 3 at first 3-MBS and then 4-MBS, etc. occur 
as the lowest excitations. In the ID case our DMRG data for the 3-magnon (octupolar) phase can be described with 
high precision by an expansion up to the second order around the critical point a 3 

h 3 s M (a) = h 03 + h' 3 {a - a 3 ) + 0.5h' 3 '(a - a 3 ) 2 + (S2) 

where h 03 = 0.1009384, h' 3 = 1.1511, and h% = 2.95699. Eq. (S2) is valid for a 3 > a > a 4 w 0.284533. A similar 
simple but semi-analytic solution, only, can be derived also for the 4-magnon phases based on a fit of our DMRG 
results. 

The inspection of Fig. S3 shows that in the limiting case a = 1 at sizable IC, say for J- lc / ' J2 > 1/3, the saturation 
field H s is also dominated by the IC. Note that a similar effect for a sizable IC appears also for a > 1. But it becomes 
weaker with increasing a. 



THE 3D SPIN-HAMILTONIAN 

The spin- 1/2 Heisenbcrg Hamiltonian reads 



R 



Ji ^Sr • Sr+i,) + J2 (Sr • £>R + 2b) + J\c ^ (Sr • Sr +i .^ + J' 1C ^2 • S R+r '^ — /i ^H • Sr^ 



(S3) 
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FIG. SI: (Color) Left: General phase diagram for frustrated chains with an AFM IC fi% of the diagonal type as realized in 
L12CUO2 (see Figs. 1 and 3 of the main text) but for arbitrary reciprocal frustration parameter 1/a = | J\\j Ji- The symbol ■ 
denotes our DMRG results and the symbol o stands for the exact hard-core boson approach. Notice a second incommensurate 
phase INC2. Its border line to the commensurate (C) phase approaches asymptotically the value of 1/a = 1 as shown by 
the dashed vertical line. For the behavior along the basis line, i.e. the ID limit, see Fig. S2 and Refs. 15 and 16 in the 
main text. The line separating the two incommensurate one-magnon phases INC1 and INC2 represents a first-order transition. 
Right: Enlarged contour plot of the dimensionless saturation field h s = g^i-B.H s /\ J\ | vs. the in-chain frustration rate a and the 
interchain interaction /J2 in the vicinity of parameters for Li2Cu02- The boundary between the INC1 and the 3-magnon phases 
is obtained by our DMRG calculations. Notice the slightly non-equidistant height differences in the contour lines of constant 
h s in the INC1 region as compared to equidistant ones in the C-region as well as the changed axis notation as compared with 
the left part. 



where R is the position vector of a Cu site in the lattice, r = (3b ± a ± c)/2, r' = (b ± a ± c)/2 a, b, c are lattice 
vectors. 

The one- and two-magnon excitation spectrum for magnetic fields H > H s can be found exactly at T = 0. The 
classical fully spin-polarized ground state is also an exact eigenstate of the quantum Hamiltonian. The spin-deviation 
(magnon) operators for the spin-iy2 Heisenberg Hamiltonian are Pauli (hard-core boson) operators, in terms of which 
the Hamiltonian is rewritten as 
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(S4) 
(S5) 

(S6) 



where [i = g/j, B and the z-axis is parallel to H. The transverse part of the Heisenberg Hamiltonian H (|S3[) defines 
the one-particle hoppings in Hq (|S5[) . the Ising part contributes the interaction (|S6[) and the on-site energy value u>o- 
Note that a negative (ferromagnetic) exchange in (|S6[) acts as an attraction between the magnons. In ID J\- J2-chains 
with Ji < and J2 > this leads always to the formation of multimagnon bound-states. In 2D and 3D lattices 
the bound-states formation is possible only if the attraction overcomes the kinetic energy. Besides on the strength of 
the AFM interchain interaction it depends also on the number of available interchain neighbors which in the case of 
diagonal IC is doubled as compared with the frequently studied but much simpler case of perpendicular IC of unshifted 
chains (see Figs. 1 of the main text and of Ref. 17 therein). Therefore, the case of a 3D chain arrangement as in 
Li2Cu02 is especially effective in deconfining the multimagnon bound states and suppressing a spiral type 1-magnon 
state as the phase INC1 and INC2 considered above. 

The Hamiltonian (|S4[) conserves the number of magnons. So, the excitation spectrum may be found separately for 
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FIG. S2: (Color) Saturation field h s vs. frustration ratio a of a single chain according to our DMRG. Notice the change in the 
slope indicated by the vertical dashed lines. The octupolar (OP or 3-MBS) is realized for 0.2845 < a < 0.3676 in between the 
hexadecupolar (HP or 4-MBS) and the quadrupolar (QP or 2-MBS) regions. Dashed line: DMRG results (each dash contains 
about 5 data points). Red solid line: Analytical curve from the exact 2-magnon solution obtained from the hard-core boson 
representation of the general spin Hamiltonian (valid exactly in the one- and two-magnon regions) given in Eq. (SI). The 
question mark in the most left yellow region close to the critical point a = 1/4 stands for numerically still unresolved but 
expected higher order multipolar phases (see text). 



one-, two-, three- etc. magnon sectors of the Hilbert space. The saturation field is determined by the condition of the 
stability of the ground state. It becomes unstable when the frequency of a certain excitation vanishes. 



THE TWO INCOMMENSURATE ONE-MAGNON PHASES INC1 AND INC2 

The energy of one-particle (1-magnon) excitations is particularly simple. For the isotropic model and J( c = it 
reads 



- /-i B H s + Ji (cos q b - 1) + J 2 (cos 2g& - 1) + 4 J ic ( cr cos cos ^ cos - 1 



(S7) 



where we have retained the AFM Brillouin zone, a = ±1 enumerates the two branches of the spectrum. The extrema 
of the one-particle spectrum are given by the equation Vu q = 0, which gives for our model q a = q c = 0, and 



. Qb 
sin — 
2 



„3 lb 



Ji cos — + 4 J 2 2 cos J — — cos — J + 3cr J ic 4 cos — - 1 



,2 lb 



= 0. 



(S8) 



The extrema for the two branches occur at = and at incommensurate values of q b = qt.a, i = 1,2, 3, where the 
expression in square brackets vanishes 
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In the dipolar phase, the saturation field H s is determined by the condition that the energy of the lowest one-magnon 
excitation w q . m i n vanishes. Then, from Eq. (|S7[) the dimensionless saturation field h s = gfj, B H s /\J\ \ is determined by 



the equation 



(cosqv 



1) — a (cos 2q 
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(S13) 



For an isolated chain (fa = 0) the two branches are degenerated and both reach the minimum at go(0) = arccos(l/4a). 
For finite J; c and a < 0.57 (the INC1 phase) the minimum occurs at of the branch a = —1. With the increase 

of fa the minimum shifts towards the centrum of BZ, which is reached at fa = /? 2 r (a) = ~ l)/9- 

Now we note that < fa < (3% < 4/9 holds for all frustration values 1/4 < a < 0.57. Then the expression (|S13|) 
may be expanded in powers of /?" — /3 2 
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(S14) 

(S15) 
(S16) 
(S17) 



For fa > (3" the minimum persists at qb,-i = 0, and the saturation field is simply given by h s — N- lc fa = 8fa. For 
Li 2 Cu0 2 one has fa w 3/76 w 0.03947 > /3f = 0.0364. For 0.57 < a < 1 the minimum occurs at q 1>+1 of the 
branch a = +1 for j3 smaller a certain value /3 2 r (a) > (4a — l)/9 where it jumps at qb.-i = 0. So, for this range of 
the frustration ratio a, the transition between the commensurate C-phasc with FM in-chain ordering and the new 
incommensurate INC2 phase is of first order. The corresponding phase boundary may be found from the condition 
w+i(g r i i -|-i) = w_i(0), and it is given by the following parametric function 



a(t) 



l-t + 



2(t-2)(2t-l) 



At 



t(t-l) 



2(2t-l) 2 
2i-7 



4ta(t) (2t - 1) + 1 - 2 
(2t — 7) (2i — 1) 



(S18) 



(S19) 



Finally, for a > 1 the minimum resides at (71,4-1 for all values of IC. The resulting phase diagram is shown schematically 
in Fig. SI (left). 



EXACT THERMODYNAMIC CONSTRAINTS INVOLVING THE SATURATION FIELD 

To demonstrate the general applicability of our proposed efficient thermodynamic method to extract exchange 
integrals from simple and "cheap" measurements, we present a collection of simple but rigorous constraints valid for 
situations where to first approximation a chain-like compound can be described by the adopted isotropic J x - J 2 -model 
for a single chain with one main interchain coupling, only, but with various geometries of the chain arrangement. 
Generalizations to several interchain couplings and to models with spin anisotropy are straightforward and will be 
given elsewhere. 

With the experimental Curie- Weiss temperature C w ( that governs the high-T spin susceptibility x(T) ~ 1/(T — 
©cw)) which is simply related to the Hamiltonian parameters by 

e cw w -0.5 [Ji + j 2 + o.5iv ic ( j ic + j/.)] (S20) 

the whole set of interchain couplings between NN chains, including in particular J lc as well as all (weak) additional 
terms such as J[ c etc. can be replaced by Q C w and Eq. (1) of the main text. As a results the we find a useful empirical 
constraint for the in-chain exchange couplings J\ and J 2 . 

. 9M Bff s +4e cw = 2|Ji|(l-a). (S21) 



b' 



Using the experimentally measured saturation field for Li2Cu02 of about 55.4 T as reported in the main text, we 
predict a ferromagnetic Curie- Weiss temperature of about 55 K to be derived from high temperature (T > 400 K) 
susceptibility measurements not yet performed. For a strictly perpendicular NN IC, only, (see e.g. Fig. 1 in Rcf. 17 ) 
the strong influence of the interchain coupling on the in-chain ordering is much reduced and as a result the 3D (2D) 
C-phase is missing, i.e. there is only one critical IC separating the single INC phase from the MP phases. In particular, 
the two-magnon phase shows a round maximum 0,01 in the 3D case near l/a = 0.906 instead of the cusp seen in 
figure SI (left) near l/a = 1.754. In the neighboring single INC-phase h s reads 

/i s = 2a-l + 0.125/a + 7V ic J ic , (S22) 

where N[ c =2 (4) for a 2D (3D) chain arrangement, respectively. Then Eq. (|S21[) is replaced by 

g^H s + 49 cw = | Ji| [1 + 0.125/a] . (S23) 

In the two magnon-phase for weak IC one has 

g f i B H s + 29 cw = | Jul [a + 0.5/(1 + a)} (S24) 

(note the new prefactor of 2 in front of 9 C w)- Similar simple constraints can be derived for any type of interchain 
couplings. They provide another part for the proposed thermodynamic approach to extract microscopic exchange inte- 
grals from few experimental numbers, only, instead to fit for instance approximate expressions for complex quantities 
such as the sussceptibility or the magmetic spefic heat in more or less extended wide temperature intervals. 

SOME DETAILED DMRG AND HARD-CORE BOSON RESULTS 

Finally, we present numerical DMRG data to demonstrate the extremely fast convergency with respect to the 
number of chains taken into account. Notice that here also a weak NN interchain coupling /3i has been taken into 
account. The DMRG result is very close to the analytical result of 0.350877192982456 given by Eq. (1) in the main 
text. 

TABLE SI: Saturation field h s at Ji = -1, a = 0.332, fa = 3/76, fa = 1/228. fa and fa are multiplied by four and two for 
2-chain and 4-chain clusters, respectively. 



L 




single chain 


2-chain cluster 


4-chain cluster 


8-chain cluster 


16 





.0610480058942 





,350877192884 





350877192884 


0.350877192884 


20 





.0618650864275 





,350877192884 





,350877192884 


0.350877192884 


24 





.0616475064797 





,350877192884 


0. 


,350877192884 


0.350877192884 


28 





.0617013393558 





,350877192884 





,350877192884 


0.350877192884 


32 





.0616885713591 





,350877192884 





,350877192884 


0.350877192884 


36 





.0616916345343 





,350877192884 


0. 


,350877192884 


0.350877192884 


40 





.0616909145148 





,350877192884 





,350877192884 


0.350877192884 


44 





.0616910784629 





,350877192884 


0. 


,350877192884 


0.350877192884 


48 





.0616910423378 





,350877192884 


0. 


,350877192884 


0.350877192884 


96 





.0616910487270 





,350877192884 








144 


.0616910487247 





,350877192884 









In order to illustrate another aspect of the accuracy of our DMRG calculations, we briefly discuss the difficult gap 
problem at H = of a single frustrated chain with a FM J\ raised by recent field theory approaches 0-0] • These 
theories predict a gap A of unknown magnitude. Since at present there are no hints from numerical studies for such 
a behavior, it is generally accepted that such a gap should be very small [§]. According to our DMRG calculations 
we estimate 10 5 1 J\\ as an upper bound. For the case of Li 2 Cu02 under consideration this corresponds to 0.01 to 
0.02 K, only. Anyhow, such a tiny gap, if it at all exists, would be of academic interest, only. 

Now we present some numerical results from exact relations obtained in the hard-core boson approach. In Fig. S3 
several examples for the IC dependence of the saturation field h s are presented. Similar dependencies are obtained for 
arbitrarily large in-chain frustration ratios a. In the 3D case with a finite antiferromagnetic J; c the saturation field 
h s is always affected by the former in contrast to the in-chain coupling which can drop out exactly what happens in 
the commensurate C-phasc. 
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FIG. S3: (Color) Saturation field for different in-chain frustration rates a = 0.5, 0.75, and 1.0 vs. IC /?2/q = J\c/ J2- The results 
of the exact hard-core boson approach are shown by symbols; the red solid line shows the one-magnon result, which becomes 
exact for 02 > f$£'~ = J^' _ / 1 -7i I as explained in the main text (see also Figs. 3-5, therein). The green dashed line shows the 
value, which is reached in the commensurate C-phase, where the dependence on Ji, J2 vanishes identically; the magenta thin 
line shows an approximation to the exact solution in the 2-magnon phase. 
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